Time-dependent approach to many-particle tunneling in one-dimension 
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Employing the time-dependent approach, we investigate a quantum tunneling decay of many- 
particle systems. We apply it to a one-dimensional three-body problem with a heavy core nucleus 
and two valence protons. We calculate the decay width for two-proton emission from the survival 
probability, which well obeys the exponential decay-law after a sufficient time. The effect of the 
correlation between the two emitted protons is also studied by observing the time evolution of the 
two-particle density distribution. It is shown that the pairing correlation significantly enhances the 
probability for the simultaneous diproton decay. 
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I. INTRODUCTION 

The quantum tunneling of a system with intrinsic de- 
grees of freedom, or of manyparticles, is an important 
subject of modern physics [I-Q. In nuclear physics, typ- 
ical examples include heavy-ion fusion reactions at sub- 
barrier energies^, Bi spontaneous fissionj^], alpha and 
heavy-cluster decaysflj, and stellar nucleosynthesis [ll|, 
[l2| . In heavy- ion fusion reactions, for instance, it has 
been well recognized that the couplings of the relative 
motion between the colliding nuclei to several collec- 
tive motions enhance the tunneling probability of the 
Coulomb barrier, thus increasing the fusion cross sec- 
tions, by several orders of magnitude as compared to a 
prediction of a simple potential model 0, Q • Neverthe- 
less, it has still been a challenging problem to under- 
stand the many-particle tunneling from a fully micro- 
scopic view. For instance, even though there have been 
several attempts [Tsl - flU , alpha decays have not fully 
been understood microscopically with sufficient accuracy. 

Recently, two-proton {2p) radioactivities have been ex- 
perimentally observed for a few proton-rich nuclei outside 
the proton drip-line, such as "^^Fe [Te'-'Ts^ and ^Be [UMI, 
and have attracted much attention [22- 25] in connection 
to e.g., the dinucleon correlations |26l - |29l |. This is a phe- 
nomenon of spontaneous emission of two valence protons 
from the parent proton-rich nuclei, in which the emission 
of one proton is energetically forbidden. Notice that an 
analogous process of the two-proton radioactivity, that is, 
a two-neutron decay has also been observed recently for 
^^Be'sO]. These two-nucleon emission decays may pro- 
vide a useful testing ground of many-particle tunneling 
theories. 

A primary task for a many-particle tunneling decay 
is to investigate the effect of interaction or correlation 
among emitted particles on decay properties such as the 
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decay width and the survival probability. In the case 
of 2p decay, this corresponds to the pairing correlation 
between the valence protons. Because of the strong pair- 
ing correlation in proton-rich nuclei, 2p emitters are ex- 
pected to have an even number of protons outside the 
proton drip-line. Incidentally, it has been well known 
that the pairing correlation plays an important role in 
two-neutron transfer reactions |3ll433l |. 

The quantum tunneling decay phenomena can be stud- 
ied either with the time- independent approach (lol[2^[25j 
or the time-dependent approach j34r 39] . In the time- 
independent approach, one seeks e.g., a Gamow state, 
which is a purely outgoing wave outside the barrier. The 
imaginary part of the energy of the Gamow state is re- 
lated to the decay width, while the real part corresponds 
to the resonance energy. On the other hand, in the time- 
dependent approach, one first modifies the potential bar- 
rier so that the initial state can be prepared as a bound 
state of a confining potential. The confining potential is 
then suddenly changed to the original barrier, and the 
initial state evolves in time. The decay width can be ob- 
tained from the survival probability of the initial state. 
An advantage of the time-independent approach is that 
the decay width can be calculated with high accuracy 
even when the decay width is extremely small (40j . An 
advantage of the time-dependent approach, on the other 
hand, is that it provides an intuitive way to understand 
the tunneling decay, even though it may be difficult to 
apply it to a situation with an extremely small decay 
width. This approach may provide a useful means to ex- 
plore the mechanism of many-particle tunneling decay, 
though it has so far been applied only to two-body de- 
cay^ph enomena, such as a decays and one-proton decays 

In this paper, we extend the time-dependent approach 
to two-proton emissions and discuss the dynamics of a 
two-particle tunneling decay. To this end, we employ the 
one-dimensional three-body model [4lj , which consists of 
a heavy core nucleus and two valence protons. We solve 
the time-dependent Schrodinger equation by expanding 
the wave function on a basis with time-dependent expan- 
sion coefficients. The decay width is then defined from 
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the survival probability. We shall study the time evo- 
lution of the survival probability as well as the density 
distribution. We shall also discuss the role of pairing 
correlation in the two-proton decay. 

The paper is organized as follows. In Sec. II we de- 
tail our formalism for the time-dependent approach to 
the two-proton decays. We show the results of the cal- 
culations in Sec. III. It will be shown that the decay 
width converges to a constant value after a sufficient time 
evolution, indicating that the decay rate follows the ex- 
ponential law. With this method, the time-evolution of 
a quasi-stationary 2p state, namely the density and the 
flux distributions, can be visualized. Using them, we 
shall discuss an important role of the pairing interaction 
between the two protons in the decay process. Finally 
we summarize the paper in Sec. IV. 
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II. FORMALISM 

A. One-dimensional three-body model 

We consider a one-dimensional three-body system with 
two valence protons and the core-nucleus whose atomic 
and mass numbers are Zc and Ac, respectively. Neglect- 
ing the recoil kinetic ener gy o f the core nucleus, the three- 
body Hamiltonian reads |41| . 



H = h{xi) + h{x2) +Vpp{xi,X2), 
2m dx"^ 



V{x), 



(1) 
(2) 



where m is the nucleon mass, xi and X2 are the coordi- 
nates of the valence protons with respect to the core nu- 
cleus. V{x) is the potential between a valence proton and 
the core, whereas Wpp(a;i,X2) is the interaction between 
the two valence protons. The core-proton potential. 



V{x) = Kuci(2;) + 14oui(a;), 



(3) 



consists of the nuclear part V^uci and the Coulomb part 
Vl:oui- For the nuclear part, we take a Woods-Saxon form. 



Kucl(2:) = - 



Vn 



1 -)_ f,i\x\-R)/a ■ 



(4) 



For the Coulomb part, we employ a soft-core Coulomb 
potentialjU-lill, that is. 



Vcoul(2:) 



\J\P- + x^ 



(5) 



Fig. 1: (Color online) The core-proton potential used in our 
calculations. The initial state at t = is constructed with a 
modified potential Vmod,(x) given by Eq. ((7|, which is shown 
by the dashed line. For t > 0, the potential is changed to 
the original potential, Vix), given by Eq. Q, as shown by 
the solid line. We also show the wave function for the bound 
state of the modified potential at e = 2.71 MeV. 



The interaction i'pp(a;i,X2) induces the pairing corre- 
lation between the two valence protons. In this work, 
we adopt a density-dependent contact interaction of the 
surface type[4l], that is. 



{xx,X2) = 1 - 



1 



1 + e(|2|--R)/a 



i(xx-X2), (6) 



where g is the strength of the interaction and x = 
(x\ 4- X2)/2. The density dependence is introduced with 
the Woods-Saxon form [R and a are the same as those in 
Eq.©). Notice that in the limit of — >■ oo, this inter- 
action becomes a pure contact interaction, ^g5{xi — X2). 
For simplicity, we neglect the Coulomb interaction be- 
tween the two protons. We have confirmed that, as long 
as the one-dimensional three-body system is concerned, 
its effect on the decay properties can be well taken into 
account by somewhat reducing the strength g (see also 



Refs. HEllil). 



It is important to notice that a one-dimensional delta 
function potential v{x) ~ ~gS{x) always holds a bound 
state at Epp = —mg'^ /Ah? for a two-proton system even 
with an infinitesimally small attraction g |47| . This is in 
contrast to a three-dimensional system, in which a bound 
state exists only with a strong strength g of an attractive 
contact interaction. 



In this paper we take Vb = 46.5 MeV. The radius R and 
the surface diffuseness parameter a in the Woods-Saxon 

1/3 

potential, Eq. are taken to be i? = 1.27 Ac fm and 
a = 0.67 fm, respectively. We arbitrary take Ac — 60 
and Zc — 30 for the mass- and atomic- numbers of the 
core nucleus, while we use 6 = 2.0 fm in the Coulomb 
interaction. 



B. Time-dependent method 

In order to describe the 2p tunneling, we employ the 
time-dependent method. The first step is to prepare the 
initial 2p state which is confined inside the potential bar- 
rier, as in the two-potential method developed by Gurvitz 
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et al. [48l - l50j . To this end, we modify the core-proton 
potential V{x) to a confining potential Vj^odix) defined 



Vtnod(a;) 



Vix) \x\ < \Xc\ 

V{Xc) \x\ > \Xc\. 



(7) 



In Fig. [U we show the confining potential, Vmod, and the 
original potential, V, together with the wave function 
for the bound state of Knod- The position Xc can be 
chosen arbitrarily as long as V{xc) is larger than the 
resonance energy, although the accuracy will be improved 
if Xc is chosen so that V{xc) is as close as possible to the 
resonance energy [sOl- When this condition is satisfied, 
the modified potential Vmod holds a bound state, which 
resembles the resonance state of the original potential. In 
this paper, we choose Xc = 10 fm, which yields the bound 
state at e=2.71 MeV for the modified potential. The 
value of the modified potential at x = Xc is Vynodixc) = 
4.21 MeV, whereas the barrier height is Vb — V{x = 
8.1 fm) =4.75 MeV. 

We solve the single-particle (s.p.) states of the modi- 
fied Hamiltonian, 



2m dx'^ 



Knod(a;), 



(8) 



(9) 



In this paper, we assume that all the s.p. states with 
negative energy, that is, e„ < 0, are occupied by the core 
nucleus. Therefore, there is only one bound state in this 
potential at e = 2.71 MeV. The other positive energy 
states are in the continuum spectra, which we discretize 
with a box of Xbox = ±120 fm. We have confirmed 
that the bound state at e=2.71 MeV corresponds to a 
resonance state of the original potential, whose energy is 
stabilized against a variation of Xbox |5ll - [53| . 

Using these s.p. wave functions, one can obtain the 
eigenfunctions of the modified three-body Hamiltonian, 

Hmod{xi,X2) ^ Kaod{xi)+hayod{x2)+Vpp{xi,X2)., (10) 

as 

*fc(a;i,X2)= ^ aW^^$„,„,(xi,a;2), (11) 

rii<Ti2 

where we exclude those s.p. states occupied by the core 
nucleus. Here, 3>nin2(2;i,a;2) is defined as 
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y/2(1 + 5^,112) 

X \S = Q). (12) 

Because we use spin-independent interactions in the 
Hamiltonian, the total spin S of the two protons is a 



good quantum number. We set it to be zero (that is, the 
spin-singlet state) . The spatial part of the 2p wave func- 
tion is therefore symmetric under the exchange of xi and 
X2 ■ The coefficients a^n^n^ in Eq. ([TTjl are determined by 
diagonalizing the Hamiltonian matrix for the modified 
Hamiltonian, H^^od- 

The next step is to carry out the time evolution for t > 
starting from the lowest eigenfunction of the modified 
Hamiltonian, i.e., ^k=o in Eq. (jlip at the eigenenergy 
Eq. That is, 



^{t = 0,a;i,2;2) = ^o{xi,X2) 



(13) 



Excluding the core-occupied states in the expansion, we 
have confirmed that there is only one bound state for 
Hm^od in the energy region between and 2V-aiod{xc)- 
Similarly to a s.p. resonance state, the lowest state cor- 
responds to the three-body resonance state of the original 
Hamiltonian H . The energy corresponds to the reso- 
nance energy of the 3-body system, that is, the Q-value 
for the 2p decay. 

We then solve the time-dependent Schodinger equation 
with the original Hamiltonian H , 



, d ^ , , 

ih— Wit, Xi,X2) 

at 



H^{t,xi,X2), (14) 
{H„,od + AV)^{t,xuX2), (15) 



where 



AV = V{xi) + V{X2) - Knod(xi) - K.od(a:2), (16) 

is the difference between the original and the modified 
potentials. We expand the time-dependent 2p wave func- 
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Fig. 2: (Color online) Overlaps between the initial state, 'Po, 
and the eigenfunctions 'i'k of the original three-body Hamilto- 
nian as a function of the corresponding eigenenergies, Ek . The 
solid and the dashed lines correspond to the cases of (; = 20 
and MeV • fm, respectively. Notice that the initial state is 
the lowest eigenstate of the modified Hamiltonian, Hmod, with 
the eigenenergy of 4.56 MeV (5.42 MeV) for = 20 MeV • fm 
(g=0 MeV-fm). 
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tion with the eigenfunctions of the modified Hamiltonian, 
that is, ^fc(xi,a;2) given in Ea.([TT|). as 



with the initial condition of 

Ck{t = 0) = 6k,o- 



(17) 



(18) 



Substituting Eq. ([T7|) into ([T5|) and using the orthogo- 
nahty of ^'fe, we obtain the differential equation for the 
expansion coefficients Ci{t), 



ih 



dci{t) 
dt 



{■^^\H\ *) (19) 

''kit) {EA.k + I I *fc)) . (20) 



Using the wave function so obtained, one can compute 
the survival probability, Ps{t), and the decay width, F, 

as 



PS) ^ \{^o\^m' = \co{i)\\ 

r . 



(21) 
(22) 



When the survival probability is an exponential function 
of i, the decay width F becomes a constant. In the next 
section, we will show that Ps (t) indeed has an exponential 
form after a sufficient time evolution. 



III. RESULTS 

A. Decay energy and width 

Before we numerically solve the time-dependent 
Schrodinger equation, let us first investigate the overlaps 
between the initial wave function and the eigenfunc- 
tions of the original Hamiltonian H . That is. 



Qk = \{^k\^o)\', 



(23) 



where \'^k) is the eigenfunctions of the original Hamilto- 
nian satisfying 



(24) 



Figure [2] shows the overlaps Qk as a function of Ek for 
g=0 and 20 MeV-fm. The initial state is fragmented 
over several eigenfunctions of the original Hamiltonian, 
and thus forms a wave packet which evolves in time. 
As one can see, the fragmentation of the initial state 
is small, where the energy spreading corresponds to the 
decay width. 

Let us now numerically solve the time-dependent 
Schrodinger equation. To this end, we use a time mesh 
of At = 0.01 fm/c. Figl3] shows the survival probability 



CO 
O 



CO 
> 

U3 



1 

0.9 
0.8 
0.7 

0.6 h 
0.5 



> 

CD 



10 
8 
6 
4 
2 



1 1 1 1 1 1 r 

(a) g=0(MeVfm) 

g=4(MeVfm) 

g=18 (MeVfm) 

9=32 (MeV fm) 



(b) 



9=0 (MeV fm) 
9=4 (MeV fm) 
9=18 (MeV fm) 
9=32 (MeV fm) 



" /'' 

I 



200 400 600 800 1000 12001400 
t (fm/c) 



Fig. 3: (Color online) (a) The survival probability as a func- 
tion of time t defined by Eq. (|2H) . (b) The decay width defined 
by Ea. (|22|l . These are plotted for several values of g indicated 
in the figure. 



and the decay width defined as Eqs. (PT|) and as a 
function of time t for several values of g. One can see 
that the decay width converges to a constant value after 
sufficient time-evolution, that indicates the exponential 
decay- law, Ps{t) — g^*'"*/''. Notice that the converged 
values for the decay width with this model Hamiltonian 
are in the same order as the experimental width for ^Be 
and ^''Ne [l9l - [2ll . Isl] . At shorter period, the decay width 
shows a transient behaviour (34j38| . That is, the sur- 
vival probability behaves like a parabolic function of i, 
whereas the decay width increases linearly [s^ . 

The dependence of the decay width on the strength 
of the pairing interaction is shown in Fig. 2] (b). We 
also show in Fig. |4] (a) the decay energy Eq (that is, 
the eigenenergy of the modified Hamiltonian given by 
Eq. (|10p ) and the asymptotic kinetic energy of diproton, 
i?rci, defined as 



E-cc\ — Eo + B-p 



(25) 



where _Bpp — mg'^ /Afi? is the binding energy of a dipro- 
ton. The decay width is estimated at i = 1200 fm/c, 
where it has been well converged (see Fig. [S]). 

The decay width F first decreases as a function of 5, 
despite that the diproton kinetic energy E-c^i increases. 
This should be related to the decrease of the decay energy 
-Bo, indicating that the sequential two-proton emissions 
is the main decay mechanism in this region of g even 
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Fig. 4: (Color online) (a) The decay energy Eo of the three- 
body system as a function of the strength of the pairing in- 
teraction, g. The dashed line indicates the asymptotic kinetic 
energy of a bound diproton. (b) The decay width estimated 
at t = 1200 fm/c. 
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(b) t=300 (fm/c) 
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though two protons arc bound in this one-dimensional 
model. For 5 > 18 MeV-fm, on the other hand, the decay 
width increases. This is consistent with the increase of 
-Brei, suggesting that the direct diproton decay, that is 
the emission of a deeply bound diproton, is the main 
mechanism in this region. 

The transition from the sequential to the diproton de- 
cays will be clarified more in the next subsection. 



Two-particle density and flux distributions 
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In order to confirm a transition from a sequential to a 
simultaneous decays discussed in the previous subsection, 
we next discuss the time evolution of two-particle density 
distribution, 

p(<,Xi,X2) =1 *(t,Xi,X2) P . (26) 

We also analyze the flux distribution defined as, 

_(^*.___*) ,,.1,2). 

(27) 

Note that the 2p density is normalized as 

/oo 
p{t,xi,X2)dxidx2 = 1- (28) 
-00 

Figs 15] and [HI show the two-particle density and flux dis- 
tributions, respectively, for g = 20 MeV • fm at i = 0, 300 



Fig. 5: (Color online) The time-evolution of the density 
distribution p{t, xi,X2) in the two-dimensional (3:1,2:2) plane 
calculated with the pairing strength of g—20 MeV ■ fm. The 
panels (a), (b), and (c) correspond to the density att = 0, 300, 
and 600 fm/c, respectively. 



and 600 fm/c. The corresponding quantities for .g = 
are also shown in Figs 17] and 15] The flux distributions 
are plotted in a form of vector at each value of {xi,X2) 
in the two-dimensional {xi,X2) plane, where the core nu- 
cleus is located at the origin. Note that there is no flux 
distribution at i = because the initial wave function 
^'(t = 0) can be taken to be real, and we do not show it 
in the figures. 

Figs. El^a) and 17][a) show the density distribution for 
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Fig. 6: (Color online) The flux distributions j{t, x-i,X2) = 
ji{t,xi,X2)e.i + j2{t, xi, X2)e2, where is the unit vector in 
the two-dimensional {xi,X2) plane. These are obtained with 
g=20 MeV-fm at 300 fm/c (Fig. 6(a)) and 600 fm/c (Fig. 
6(b)), and are plotted in arbitrary units. 



the initial 2p state, which is confined within the modified 
potential, Knod- Because of the pairing correlation, the 
initial density for g = 20 MeV • fm has an asymmetric 
form, with the peaks along xi = X2 being higher than 
those along xi — —X2^^- The peaks along the xi = x^ 
line, that is, in the first and third quadrants of these 
panels, correspond to a compact diproton cluster. On 
the other hand, the peaks along the x\ — ~X2 hne {i.e., 
in the second and the fourth quadrants) correspond to a 
configuration in which two protons are located opposite 
to the core nucleus. If we discard the pairing interac- 
tion, the density distribution has four symmetric peaks, 
as shown in Fig. [3 that is, the probability in the first 
and third quadrants is the same as that in the second 
and fourth quadrants [41]. 

The effect of pairing correlation is apparent also during 
the time evolution. In the presence of the pairing corre- 
lation, the extension of the two-particle density along the 
xi — X2 hne increases significantly, although the exten- 
sion along the xi = or X2 = lines is not negligible. 
This is in marked contrast with the uncorrelated case 
shown in Fig. [71 in which the two-particle density ex- 
pands democratically. That is, in the uncorrelated case. 
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Fig. 7: (Color online) Same as FigO but for g=0 MeV • fm. 



the probability of emission of the two protons in oppo- 
site directions is equal to that in the same direction. The 
flux distribution shown in Figs. [6] and [8] also indicate the 
same behaviour. 

In order to investigate the time evolution more quan- 
titatively, we divide the {xi,X2) plane into four regions 
shown in Fig. [9l That is, (i) the region of xi > 16 fm 
and X2 > 16 fm, as well as the region of xi < — 16 fm and 
X2 < — 16 fm, (ii) the region of xi > 16 fm and a;2 < — 16 
fm, as well as the region of xi < —16 fm and X2 > 16 fm, 
(iii) the region of — 16 < Xi < 16 fm and \x2\ > 16 fm, 
as well as the region of —16 < X2 < 16 fm and \xi \ > 16 
fm, and (iv) the rest in the {xi, X2) plane. At each time. 
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Fig. 8: (Color online) Same as Fig[6l but for g—0 MeV ■ fm. 



we integrate the two-particle density distribution in each 
region, 



Pkit) 



p{t,xi,X2)dxidx2, {k — l^A). (29) 



region k 



The time evolution of these partial probabilities is shown 
in FigHnifor g = 0, 20 and 32 MeV - fm. We also show the 
total decay probability, that is, 1 — P4 = Pi + P2 + ^3- 

We first discuss the behaviour for the uncorrelated case 
shown in Fig. [10] (a). In this case, the dominant process is 
the decay into the third region, P3 , which corresponds to 
an emission of one of the valence protons while the other 
proton remains inside the core-proton potential. As there 
is no active s.p. bound state in the present core-proton 
potential, the second proton mainly occupies the s.p. res- 
onance state. This resonance state eventually decays and 
the second proton is emitted outside the potential after 
a sufficient time-evolution. This is nothing but the se- 
quential two-proton decay, and can be clearly seen in the 
probabilities Pi and P2, that exist only at t >; 600 fm/c. 
Notice that Pi and P2 are identical to each other, since 
the second proton is emitted into either the left or the 
right direction with respect to the core nucleus with an 
equal probability irrespective to the position of the first 
proton. 

For g = 20 MeV • fm shown in Fig. [TUTb). the proba- 
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Fig. 9: (Color online) The four regions in the {x-i,X2) plane 
used to calculate the partial probabilities shown in Fig. [TO] 
The boundaries of each region are at xi = ±16 fm and X2 = 
±16 fm. 



bility in the region (i) increases considerably due to the 
pairing correlation, while P3 decreases significantly. This 
partly corresponds to an emission of a bound diproton, 
that is, the simultaneous two-proton decay. Notice, how- 
ever, that P3 is still larger than Pi at i < 800 fm/c, 
and a sequential decay also coexists for this value of g. 
As we have shown in FigUl the total decay probability. 
Pi + P2 + P3, decreases compared to the uncorrelated 
case. 

When the pairing interaction is even stronger, the 
simultaneous diproton decay becomes dominant. See 
FigHni (c) for g = 32 MeV • fm. In this case, Pi is the 
dominant part of the total decay probability, except for 
the short time region, at which the high energy compo- 
nents in the initial wave function quickly escape from the 
potential barrier. The long-time behaviour in this case 
may correspond to alpha decays in realistic nuclei, for 
which a tightly bound alpha particle tunnels through the 
Coulomb barrier of the daughter nucleus. 

From these studies, it is evident that the present one- 
dimensional three-body model nicely describes a transi- 
tion from an uncorrelated case to a strongly correlated 
case for many-particle tunneling decays. 



IV. SUMMARY 

We have employed the time-dependent method and in- 
vestigated many-particle tunneling decays, particularly 
the two-proton radioactivity. To this end, we have used 
a one-dimensional three-body model which consists of a 
core nucleus and two valence protons. In order to de- 
scribe the decaying process, we first confined the two- 
proton wave function inside a confining potential. The 
confining potential was then changed to the original po- 
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time evolution of the density and flux distributions. We 
have also analyzed the partial probabilities and discussed 
the relative importance of the sequential and simultane- 
ous two-proton decays. We have shown that, for the un- 
correlated case, the sequential decay is the dominant de- 
cay process, while the simultaneous decay plays an essen- 
tial role in the case of strong pairing correlation. For an 
intermediate value of the pairing strength, we have shown 
that both the simultaneous and the sequential two-proton 
emission coexist. 

The one-dimensional three-body model which we em- 
ployed in this paper is a simple schematic model, with 
which a deep understanding of many-particle decay pro- 
cess may be obtained. One drawback, however, is that 
two protons are inevitably bound even with an infinitesi- 
mal attraction between the two protons. Even though 
our model nicely demonstrates the coexistence of the 
simultaneous and sequential decays for an intermedi- 
ate pairing interaction, in reality two protons are never 
bound in vacuum. It will be an important task to ex- 
tend the present study to realistic two-proton emitters in 
three-dimension, such as ^Be and ^^Ne nuclei. A work 
towards this direction is in progress, and we will report 
on it in a separate publication. 



Fig. 10: (Color online) The time-evolution of the partial prob- 
abilities in the regions defined in Fig. |5] for g=Q, 20, and 
32 MeV-fm. The total decay probability, Ptot = 1 - P4 = 
Pi + P2 + P3, is also shown by the dot-dashed lines. 



tential, with which the two-proton wave function evolves 
in time. We have confirmed that the survival probabil- 
ity follows the exponential decay-law after a sufficient 
time-evolution, yielding a constant decay width. We have 
found that an emission of a bound diproton is enhanced 
due to the pairing correlation, as was evidenced in the 
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